40 research outputs found
An iterative semi-implicit scheme with robust damping
An efficient, iterative semi-implicit (SI) numerical method for the time
integration of stiff wave systems is presented. Physics-based assumptions are
used to derive a convergent iterative formulation of the SI scheme which
enables the monitoring and control of the error introduced by the SI operator.
This iteration essentially turns a semi-implicit method into a fully implicit
method. Accuracy, rather than stability, determines the timestep. The scheme is
second-order accurate and shown to be equivalent to a simple preconditioning
method. We show how the diffusion operators can be handled so as to yield the
property of robust damping, i.e., dissipating the solution at all values of the
parameter \mathcal D\dt, where is a diffusion operator and \dt
the timestep. The overall scheme remains second-order accurate even if the
advection and diffusion operators do not commute. In the limit of no physical
dissipation, and for a linear test wave problem, the method is shown to be
symplectic. The method is tested on the problem of Kinetic Alfv\'en wave
mediated magnetic reconnection. A Fourier (pseudo-spectral) representation is
used. A 2-field gyrofluid model is used and an efficacious k-space SI operator
for this problem is demonstrated. CPU speed-up factors over a CFL-limited
explicit algorithm ranging from to several hundreds are obtained,
while accurately capturing the results of an explicit integration. Possible
extension of these results to a real-space (grid) discretization is discussed.Comment: Submitted to the Journal of Computational Physics. Clarifications and
caveats in response to referees, numerical demonstration of convergence rate,
generalized symplectic proo
Simulation of the Magnetothermal Instability
In many magnetized, dilute astrophysical plasmas, thermal conduction occurs
almost exclusively parallel to magnetic field lines. In this case, the usual
stability criterion for convective stability, the Schwarzschild criterion,
which depends on entropy gradients, is modified. In the magnetized long mean
free path regime, instability occurs for small wavenumbers when (dP/dz)(dln
T/dz) > 0, which we refer to as the Balbus criterion. We refer to the
convective-type instability that results as the magnetothermal instability
(MTI). We use the equations of MHD with anisotropic electron heat conduction to
numerically simulate the linear growth and nonlinear saturation of the MTI in
plane-parallel atmospheres that are unstable according to the Balbus criterion.
The linear growth rates measured from the simulations are in excellent
agreement with the weak field dispersion relation. The addition of isotropic
conduction, e.g. radiation, or strong magnetic fields can damp the growth of
the MTI and affect the nonlinear regime. The instability saturates when the
atmosphere becomes isothermal as the source of free energy is exhausted. By
maintaining a fixed temperature difference between the top and bottom
boundaries of the simulation domain, sustained convective turbulence can be
driven. MTI-stable layers introduced by isotropic conduction are used to
prevent the formation of unresolved, thermal boundary layers. We find that the
largest component of the time-averaged heat flux is due to advective motions as
opposed to the actual thermal conduction itself. Finally, we explore the
implications of this instability for a variety of astrophysical systems, such
as neutron stars, the hot intracluster medium of galaxy clusters, and the
structure of radiatively inefficient accretion flows.Comment: Accepted for publication in Astrophysics and Space Science as
proceedings of the 6th High Energy Density Laboratory Astrophysics (HEDLA)
Conferenc
Large-Eddy Simulations of Magnetohydrodynamic Turbulence in Heliophysics and Astrophysics
We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) ows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several special applications in heliophysics and astrophysics, assessing triumphs, challenges,and future directions
Recommended from our members
A Numerical Instability in an ADI Algorithm for Gyrokinetics
We explore the implementation of an Alternating Direction Implicit (ADI) algorithm for a gyrokinetic plasma problem and its resulting numerical stability properties. This algorithm, which uses a standard ADI scheme to divide the field solve from the particle distribution function advance, has previously been found to work well for certain plasma kinetic problems involving one spatial and two velocity dimensions, including collisions and an electric field. However, for the gyrokinetic problem we find a severe stability restriction on the time step. Furthermore, we find that this numerical instability limitation also affects some other algorithms, such as a partially implicit Adams-Bashforth algorithm, where the parallel motion operator v{sub {parallel}} {partial_derivative}/{partial_derivative}z is treated implicitly and the field terms are treated with an Adams-Bashforth explicit scheme. Fully explicit algorithms applied to all terms can be better at long wavelengths than these ADI or partially implicit algorithms
Recommended from our members
Comparing Linear Microinstability of the National Compact Stellarator Expriment and a Shaped Tokamak
One metric for comparing con nement properties of di erent magnetic fusion energy con gurations is the linear critical gradient of drift wave modes. The critical gradient scale length determines the ratio of the core to pedestal temperature when a plasma is limited to marginal stability in the plasma core. The gyrokinetic turbulence code GS2 was used to calculate critical temperature gradients for the linear, collisionless ion tem- perature gradient (ITG) mode in the National Compact Stellarator Experiment (NCSX) and a prototypical shaped tokamak, based on the pro les of a JET H-mode shot and the stronger shaping of ARIES-AT. While a concern was that the narrow cross section of NCSX at some toroidal locations would result in steep gradients that drive instabilities more easily, it is found that other stabilizing e ects of the stellarator con guration o set this so that the normalized critical gradients for NCSX are competitive with or even better than for the tokamak. For the adiabatic ITG mode, NCSX and the tokamak had similar critical gradients, though beyond marginal stability, NCSX had larger growth rates. However, for the kinetic ITG mode, NCSX had a higher critical gradient and lower growth rates until a/LT ≈#25; 1:5 a/LT;crit, when it surpassed the tokamak's. A discussion of the results presented with respect to a/LT vs R/LT is included
Recommended from our members
Measurements of energetic helium-3 minority distributions during ion cyclotron radio-frequency heating in the Princeton Large Torus
Ion cyclotron radio-frequency heating experiments were performed with a /sup 3/He minority ion species in a /sup 4/He majority plasma in the Princeton Large Torus. The energetic /sup 3/He ion ''tail'' was measured directly with a charge exchange neutral analyzer for the first time. Comparisons with bounce-averaged quasi-linear calculations suggest a modestly peaked radi-frequency power deposition profile. The double charge exchange process /sup 3/He/sup + +/ )plus) /sup 4/He/sup 0/ )plus) /sup 3/He/sup 0/ )plus) /sup 4/He/sup + +/ demonstrated in these measurements may be useful as part of an alpha particle diagnostic in a fusion reactor experiment. 18 refs., 4 figs